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ABSTRACT 


The  31P(d,n) 
energies  of  4.0  and 
by  the  Time-of -Flight 


S  reaction  has  been  studied  at  deuteron 
5.45  MeV.  Neutron  energies  were  measured 
technique.  The  angular  distributions  of 
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seven  levels  in  S-nucleus  have  been  measured  and  compared  with 
theoretical  predictions  of  both  the  DWBA  using  the  computer  code 
DUWK  and  the  compound  statistical  model  using  Hauser-Feshbach  code. 
Analysis  of  angular  distributions  yields  the  values  along  with 
the  absolute  spectroscopic  factors.  The  results  are  compared  with 
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shell  model  prediction  and  other  (d,n)  and  (He  ,d)  reactions. 
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CHAPTER  1 


INTRODUCTION 


The  literature  of  spectroscopic  studies  of  (d,p)  and  (  d,n) 
reactions  is  characterized  by  two  distinct  periods  of  time.  In  the  early 
period  these  reactions  have  been  analyzed  by  using  the  Plane-Wave-Born- 
Approximation  theories  (Bu51)  and  (Bh52) .  These  theories  were  able  to 
predict  spectroscopic  quantities  such  as  angular  momenta  "£"  of 
stripped  particles,  parities  "tt"  of  energy  levels  of  residual  nucleus, 
and,  in  some  cases,  a  definite  assignment  of  their  spins  "J".  This  early 
period  came  to  its  end  with  the  use  of  the  Distorted-Wave- Born-Approxima¬ 
tion  (DWBA)  in  calculating  the  angular  distributions  for  stripping  reactions. 
This  theory  is  able  to  predict  not  only  the  above  mentioned  spectroscopic 
quantities  but  also  in  most  cases  it  can  predict  correctly  the  magnitude 
of  the  main  peak  of  the  angular  distribution.  So  it  became  easy  to  extract, 
in  addition  to  the  conventional  spectroscopic  quantities  tt  and  J" , 

a  very  important  quantity,  'The  Spectroscopic  Factor  S'. 

31  32 

The  literature  of  P(d,n)  S  reaction  followed  a  track  very 
similar  to  that  outlined  above. 

31  32 

Early  investigations  of  the  P(d,n)  S  reaction  by  Calvert 

et  al  (Ca  55)  and  by  El-Bedewi  et  al.  (El  55)  have  yielded  information  about 
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the  lowest  energy  levels  in  S  .  Since  then  the  only  reported  work  on 

this  reaction  is  that  of  Ferguson  et  al.  (Fe  68) .  They  carried  out  this 

reaction  at  E,  =  4.9  MeV  and  extracted  the  Z  values  and  isobar ic  spins 

d  P 
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for  some  energy  levels  in  S.  They  did  not,  however,  report  either  the 
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absolute  magnitude  of  the  cross-section  nor  spectroscopic  factors. 

31  32  32 

Apart  from  the  P(d,n)  S  reaction,  the  S  nucleus  has 
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been  studied  by  other  reactions.  All  information  about  S  nucleus 
concerned  with  Z,  tt  ,  J  and  excitation  energies  below  9  MeV  have  been 
compiled  by  Endt  et  al.  (En  67). 

31  32 

The  above  survey  shows  that  the  absolute  magnitude  of  P(d,n)  S 
cross-reactions  and  extraction  of  spectroscopic  factors  are  missed  from  the 
literature.  So  this  reaction  has  been  done  twice  and  Z,  tt,  S,  a  (6)  has 


been  extracted  and  compared  with  those  recently  reported  by  Graue  et  al. 

3 

(Gr  68)  and  Morrison  (Mo  70)  from  (He  ,d)  reaction  at  deuteron  energies  of 
12  MeV  and  15  MeV  respectively.  The  spectroscopic  factors  are  also  compared 
to  those  obtained  by  Glaudmans  et  al.  (G1  64)  from  shell  model  calculation  in 
the  s-d  shell  region  of  nuclei. 
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CHAPTER  2 

THEORETICAL  BACKGROUND 


2 . 1  Introduction 

The  study  of  nuclear  reactions  is  of  prime  importance  in  the 
investigation  of  the  properties  of  nuclear  matter  and  of  nuclei.  A 
nuclear  reaction  results  in  general  if  a  nucleus  X  ,  say,  is  bombarded 
with  a  nuclear  particle  a  which  may  be  a  nucleus  or  a  nucleon.  Such 
a  nuclear  reaction  is  usually  written  symbolically  as: 

a  +  X  -*  Y  +  b  2.1 

or,  in  more  compact  notation  X(ayb) Y  .  This  means  that  a  nucleus  Y 
has  been  formed  by  the  bombardment  of  X  with  a  ,  and  the  outgoing 
particle  b  carries  off  the  excess  energy  originally  brought  in  by  the 
particle  a  .  In  the  particular  case  that  a  and  b  are  particles  of 
the  same  type,  Y  will  represent  the  same  nucleus  as  X  ;  it  will,  however, 
be  in  a  different  quantum  state  from  that  of  the  bombarded  nucleus  (which 
is  usually  in  the  ground  state)  for  all  cases  in  which  the  kinetic  energy 
of  the  outgoing  particle  differs  from  that  of  the  incident  particle  (in  the 
center-of-mass  system). 

Of  course  when  a  large  number  of  nuclei  X  in  a  target  are 
bombarded  by  a  beam  of  particles  of  type  a  ,  a  large  number  of  different 
nuclear  reactions  will  result.  A  given  reaction  may  be  studied  experiment¬ 
ally  by  the  detection  of  particles  b  of  a  specific  type  and  having  a 
certain  particular  energy;  the  type  of  particles  b  detected  determines 
the  nature  of  the  nucleus  Y  formed  in  the  reaction  being  studied,  and  the 
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energy  of  the  outgoing  particles  determine  the  particular  state  Y 
appropriate  to  this  reaction. 

The  information  that  can  be  obtained  from  a  study  of  nuclear 
reactions  may  be  divided  into  two  separate  categories: 

1)  Information  concerning  properties  of  nuclear  matter  in 
general. 

2)  Information  concerning  specific  properties  such  as  spin, 
parity,  spectroscopic  factors,  etc.  of  the  particular 
nuclear  state  Y  formed  in  a  given  reaction. 

As  far  as  this  work  is  concerned  we  are  interested  in  the  obtain¬ 
ing  of  information  in  category  (2);  the  use  of  nuclear  reactions  as  a  tool 
for  the  investigation  of  specific  properties  of  nuclear  states,  i.e.,  as  a 
tool  in  the  field  of  nuclear  spectroscopy.  For  this  purpose  it  is  desirable 
to  choose  the  type  of  nuclear  reaction,  or  the  conditions  under  which  a 
given  reaction  is  performed,  so  that  the  observations  made  on  such  reactions 
are  directly  interpretable  in  terms  of  the  properties  of  the  initial  and 
final  nuclear  levels  involved.  If,  as  usual,  the  properties  of  the  initial 
nucleus  are  known,  the  corresponding  properties  of  the  final  state  may  then 
be  deduced.  This  is  not  an  easy  job  in  the  absence  of  an  exact  theory  of 
nuclear  reactions  and  lack  of  enough  knowledge  about  the  actual  mechanisms 
through  which  the  reaction  of  Eq.  2.1  can  proceed. 

Efforts  in  the  theory  of  nuclear  reactionsare  largely  devoted  to 
attempts  at  finding  simple  models  that  are  capable  of  representing  one  or 
more  aspects  of  the  experimental  situation  (To  61).  Two  models  of  nuclear 
reactions  have  enjoyed  particular  success;  the  compound  nucleus  model  and 
the  direct  reaction  model.  According  to  the  compound  nucleus  model  a 


' 
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nuclear  reaction  takes  place  in  two  steps:  a)  The  target  nucleus  captures 
the  incident  particle  to  form  a  metastable  state,  the  compound  nucleus, 
b)  The  compound  nucleus  subsequently  disintegrates  to  yield  the  reaction 
products.  The  existence  of  compound  nucleus  states  of  relatively  long  life 
time  provides  an  explanation  for  the  narrow  resonances  observed  in  nuclear 
cross-sections  at  low  energies.  According  to  the  direct  reaction  model 
the  transition  from  the  incident  channel  to  the  reaction  channel  in  a 
nuclear  reaction  takes  place  in  one  step  without  the  formation  of  an  inter¬ 
mediate  state.  In  this  model  the  target  nucleus  and  incident  projectile 

-22 

system  has  a  life  time  of  the  order  of  10  sec,  assuming  an  interaction  - 

potential  depth  of  some  tens  of  MeV.  On  the  other  hand  the  compound  nucleus 

-14 

has  a  life  time  of  the  order  of  10  sec  for  an  energy  width  of  the  order 
of  a  fraction  of  an  electron  volt.  Thus  the  time  scales  of  the  two  mechan¬ 
isms  are  quite  different.  The  cross-sections  predicted  by  the  direct  reaction 
theory  do  not  show  the  sharp  resonance  characteristics  of  the  compound  nucleus 
cross-sections.  A  distinctive  feature  of  the  direct  reaction  cross-sections 
is  the  strong  dependence  on  the  direction  of  emission  of  the  reaction  products. 

In  the  following  we  will  discuss  in  some  detail  and  give  a  quick 
sketch  of  the  most  successful  theories  for  each  model  separately. 

2 . 2  The  Direct  Reactions 

The  direct  reaction  is  a  one-step  process.  This  was  first  recog¬ 
nized  by  Oppenheimer  and  Phillips  (Op  35)  in  analyzing  the  low  energy  (d,p) 
reactions  (La  35). 


. 
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Oppenheiraer  and  Phillips  explained  the  reaction  by  noting  that 
the  deuteron  is  a  loosely  bound  system  (low  binding  energy  of  2.2  MeV  and 
radius  large  compared  to  the  range  of  two  nucleon  force) ,  and  when  it 
approaches  the  target  nucleus  the  proton  is  detached  from  the  deuteron 
due  to  the  Coulomb  field,  whereas  the  neutron  is  captured.  At  low  energies 
this  is  known  as  the  Oppenheimer  and  Phillips  process 3  or  the  Coulomb 
Stripping  process 3  although  it  is  now  more  commonly  known  as  the  Stripping 
process j  at  both  low  and  high  energies.  In  the  high  energy  region  the 
(d,n)  and  (d,p)  reactions  are  equally  probable.  But  what  distinguishes 
the  direct  reaction  from  other  nuclear  reaction  mechanisms  is  that  the 
angular  distribution  of  disintegration  products  is  peaked  in  the  forward 
direction  with  a  very  small  intensity  in  the  backward  direction. 

It  was  found  that,  not  only  the  (d,n)  and  (d,p)  reactions 

but  also,  a  large  number  of  nuclear  reactions,  such  as  (°Sp)  ,  (p,d), 

3  3 

(n,d)  ,  (He  ,<*)  ,  (He  ,d)  and  many  others,  showed  the  characteristic  of 
stripping  or  the  inverse  pick-up  reactions.  Reactions  of  types  (p,p')  , 
(‘s,0'’)  and  (d,d')  can  also  be  analyzed  according  to  the  direct  process. 

In  the  use  of  nuclear  reactions  as  a  tool  in  nuclear  spectroscopy 
one  is  concerned  with  extracting  information  about  specific  levels  in  the 
final  nucleus  Y  .  The  interest  of  the  stripping  process  arises  from  the 
fact  that  when  the  absorbed  particle  forms  a  definite  quantum  state  of  the 
residual  nucleus,  it  must  take  in  a  certain  definite  orbital  angular  momen¬ 
tum  and  it  must  cause  a  definite  change  in  kinetic  energy.  These  changes 
in  orbital  angular  momentum  and  energy  are  impressed  on  the  wave  functions 
of  the  outgoing  particle,  and  a  measure  of  the  angular  distribution  of  this 
particle  can  give  a  determination  of  t  -  value  of  the  absorbed  particle  and 
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hence  the  parity  change  between  the  initial  and  final  states.  The  dependence 
of  the  angular  distribution  on  the  £-value  can  be  explained  by  the  following 
simple  classical  arguments  (Bu  57). 


Consider  a  deuteron  stripping  reaction  in  the  form 


x4 


+  d 


/+! 


+  b 


2.2 


where  d  =  b+c  ,  d  =  deuteron,  b  =  outgoing  particle,  c  =  captured  particle, 
/+1  -  xA  +  c  =  residual  nucleus  and  =  target  nucleus.  Consider  also 
Fig.  2-1,  which  represents  the  vector  relations  existing  between  the  various 

particle-momenta  immediately  after  stripping  has  taken  place.  The  momentum 

-> 

triangle  dbc  shows  how  the  original  deuteron  momentum  TiK^  is  shared 
between  the  outgoing  and  captured  particle  momenta,  TiK^  and  TiK^  respect- 


-> 


ively  in  the  laboratory  system.  Tiq^  is  the  internal  momentum  contribution 
•'  -> 

to  the  outgoing  momentum  "hK^  ,  and  vectorially  we  get: 

%  =  h-  ah 


q,,  ■ 


/k^  +  iKd  *  hh  cos  6 


2.3 


and 


->  -> 

K  = 
c 


5, 


K. 


K 


=  /  2 


K  +  h  -  2\ 


K,  cos  0 
d 


2.4 


The  angular  momentum  carried  by  the  captured  particle  into  the  target  is 

TiK  n  ,  where  n  is  the  classical  impact  parameter,  for  the  reaction  to 
c 

proceed,  TiK  n  must  be  at  least  as  large  as  ,  the  orbital  angular 

c  c 

momentum  of  the  state  into  which  the  particle  c  is  captured,  i.e., 

K  n  >  l 
c  —  c 


2.5 
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FIGURE  2-1 


Vector  Relation  For  Deuteron  Stripping 
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This  is  a  very  simple  expression  of  the  conservation  of  angular  momentum. 
Condition  (2.5)  leads  to  two  cases  of  interest: 

a)  t  =  0  ,  in  this  case  condition  (2.5)  will  be  satisfied  for  all  values 

of  ,  and  hence,  for  all  angles  of  scattering  0  .  The  main  features 

determining  the  shape  of  the  angular  distribution  will  be,  therefore,  the 
fact  that  the  deuteron  will  always  prefer  to  break-up  in  such  a  way  that 
the  vector  is  close  to  the  direction  of  motion  of  the  deuteron.  In 

this  case,  =  0  ,  the  angular  distribution  is  peaked  in  the  forward 
direction. 

b)  t  >  0  ,  in  this  case  the  condition  K  n  =  t  cannot  be  satisfied 

c  c  c 

for  small  angles  of  scattering  (since  Kc  -*  0  as  0  ->  0) .  The  scattering 
angle  0  must  increase  to  a  value  0^  such  that  K  r\  =  £  before  the 

reaction  will  proceed.  Thus  a  pronounced  maximum  is  expected  at  the  angle 
0Q  .  As  increases,  this  maximum  will  shift  out  to  larger  angles  to 

satisfy  the  condition  (2.5). 


Taking  a  more  quantum  mechanical  view  point  we  can  introduce 


restrictions  on  the  allowed  values  of  t 


If  the  initial  and  final 


nuclear  spins  are  J\  and  respectively,  and  the  captured  particle  has 


-> 


spin  ,  then 


J  =  J  .  +  l  +  S 
f  1  c  ( 


2.6 


l.  e. 


|J.  -  J  J  -  4  <  l  <  Jf  +  j,  +  i 

'i  f1  2  —  c  —  f  i  2 


2.7 


or 


<  J .  +  l 
—  1  c 


2.8 
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This  is  for  the  case  of  .  The  parity  of  a 
state  with  orbital  angular  momentum  quantum  number 
and  hence,  conservation  of  parity  requires  that 


quantum  mechanical 

l 

lc  is  (-1)  C  , 


l 

*.(-1)  C 


2.9 


where  tt  and 

l 


7T 

f 


are  the  initial  and  final  state  parities  respectively. 


Additional  nuclear-structure  information  can  be  obtained 
from  the  intensities  of  the  outgoing  particle  groups.  The  deuteron  strip¬ 
ping  reactions  are  highly  selective,  strongly  populating  only  those  states 
in  the  residual  nucleus  that  can  be  obtained  by  simply  adding  a  neutron  or  a 
proton  to  the  ground  state  of  the  target  nucleus.  Several  single  particle 
orbitals  may  be  available  for  the  captured  particle,  and  the  degree  to  which 
the  residual  nucleus  satisfies  this  condition  for  a  particular  single 
particle  orbital  (£,j)  is  called  its  spectroscopic  factor  S(£,j)  for 
that  orbital.  The  S(£,j)  is  unity  if  the  state  exhausts  the  single  particle 
strength,  as  would  occure  in  capture  by  a  closed-shell  nucleus  into  a  pure 
single  particle  state. 


A  formula  for  the  deuteron  stripping  angular  distribution  was  first 
obtained  by  S.T.  Butler  (Bu  51)  using  the  difficult  mathematical  procedure 
of  matching, in  phase  and  slope  at  the  nuclear  boundary,  the  wave  functions 
for  the  external  and  internal  regions. 


Using  alternative  methods,  several  authors  obtained  the  same  or 
similar  results  to  those  of  Butler,  e.g.  Bhatia  et  aL  (Bh  52)  used  the  Plane- 
Wave-Born-Approximation  PWBA ,  and  although  this  procedure  was  difficult  to 
justify  physically,  later  and  more  detailed  studies  (Au  53,  To  54)  showed 


< 
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that  Butler's  approach  was  in  the  final  analysis  equivalent  to  the  PWBA 
approach. 


The  PWBA  theory  was  remarkably  successful  as  a  tool  for  identi¬ 
fying  £  -values  from  measured  angular  distributions  from  light  nuclei 
and  even  gave  good  fits  to  the  shape  of  the  main  peak  in  the  angular 
distribution.  However, it  was  quickly  noted  (Fu  53)  that  the  predicted 
absolute  cross-sections  were  frequently  too  large  by  an  order  of  magnitude 
or  more.  Thus,  the  absolute  values  of  spectroscopic  factors  could  not  be 
obtained.  However,  the  failure  to  predict  absolute  cross-sections  also 
cast  some  doubt  on  the  ability  of  PWBA  theory  to  give  a  correct  account  of 
relative  cross-sections  and  their  dependence  on  energy  and  Q-value.  As 
experimental  data  for  heavier  nuclei  became  available,  where  the  Coulomb 
force  at  the  bombarding  energies  usually  available  is  too  important  to  be 
ignored,  it  became  clear  that  the  angular  distributions  also  deviated 
considerably  from  the  expectation  of  the  simple  theory.  It  is  now  known 
that  distortion  by  nuclear  scattering  and  absorption  is  important,  and  hence, 
a  plane-wave  theory  is  not  a  good  approximation. 

To  obtain  a  detailed  agreement  between  theory  and  experiment  for 
stripping  reactions  involving  light  and  intermediate  nuclei,  and  to  extract 
spectroscopic  information  in  the  heavy- elements  region,  it  was  proposed 
that  the  interaction  could  be  more  accurately  described  by  the  distorted- 
wave-Born-Approximation  (DWBA)  (To  54) .  The  DWBA  theory  takes  account  of 
the  scattering  and  absorption  of  the  incident  deuteron  before  stripping,  and 
of  outgoing  particle,  by  replacing  the  plane  waves  by  distorted  or  elastic 
scattering  waves.  In  practice,  these  are  generated  by  optical-model  poten¬ 
tials  that  reproduce  the  observed  elastic  scattering  from  the  same  nucleus 
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at  the  same  energy,  and  whose  parameters  are  therby  determined.  Consider¬ 
able  effort  has  gone  into  the  exploitation  of  this  theory,  and  numerous 
comparisons  with  experimental  data  have  been  done  and  suggested  that 
quantitatively  accurate  predictions  are  possible  (for  list  of  references 
see  e.g.  Sa  63,  G1  63  and  Le  64). 


To  account  for  the  many  features  of  the  reaction,  it  is  possible 
to  divide  it  into  three  separate  stages  (El  69) : 

a)  The  projectile  (deuteron)  moves  in  the  average  field  of  the  target. 
This  field  can  be  described  by  means  of  an  optical  potential,  and  leads 
to  a  wave  function  identical  to  the  one  describing  the  elastic  scattering 
of  deuterons. 

b)  A  nucleon  is  transfer ed  to  an  empty,  single  particle  orbital  in  the 
target.  In  may  cases,  it  is  a  good  approximation  to  consider  the  target 
core  as  unchanged  by  this  capture. 

c)  The  outgoing  particle  emerges  acted  upon  by  the  field  of  the  final 
nucleus.  This,  again,  is  an  elastic  scattering  problem. 


These  assumptions  suggest  that  the  differential  cross-section 


can  be  factorized  as: 


^  =  No  „(9)  scej) 
dco  jZ 


N a.«(0)  is  the  cross-section  for  adding  (or  removing  in  pick-up)  a 
JZ 

particle  to  a  pure  single  particle  orbital  (£,j)  >  0  is  the  direction  of 
the  outgoing  particle.  The  second  factor  S(£,j)  is  the  spectroscopic 
factor  defined  before,  which  contains  the  spectroscopic  information.  The 
single  particle  cross-section  Na  ^(0)  is  determined  by  the  parameters 
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describing  the  elastic  scattering  process. 

Let  us  now  calculate  the  angular  distribution  of  the  outgoing 
particles  in  the  reaction  X(atb) Y  .  The  Hamiltonian  may'be  separated  (Sa  65) 
into  the  Hamiltonian  of  the  colliding  pair  when  not  interacting  (well 
separated) , and  the  interaction  V  between  them, 


II 

re 

+ 

V 

+ 

T 

2.11 

oc 

OC 

OC 

=  H 

+ 

v0 

+ 

T  „ 

2.12 

3 

3 

3 

where  H  =  £Lr  +  H  and 
colliding  pair  X  and  a 
the  Schrodinger  equation 


is  the  kinetic  energy  operator  of  the 
.  The  total  wave  function  of  the  system  satisfies 


(E-H)y  =  0 


2.13 


while  the  wave  functions  of  the  internal  states  of  the  pair  are  given  by 


(E  -H  )4>  =  0 


<W*B  =  ° 


\b  \h 

Va 


2.14 


and  E  =  E  +  E  ,  and  the  total  energy  of  the  system  is  given  by: 

«  X  a 


E  =  E 

2  2 
+  iTk 

/ 

2m 

OC 

OC 

OC 

2.15 

=  E. 

2  2 
+  Ti  kf 

/ 

2m . 

3 

3 

3 

where  k  is  the  wave  number  of  the  relative  motion  of  the  X  and  a 

CC 

and  m  is  their  reduced  mass. 

oc 

The  theory  of  direct  nuclear  reactions  is  concerned  with  calcula¬ 
ting  the  transition  amplitudes  for  various  types  of  nuclear  reactions  on 
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the  basis  outlined  above.  The  direct  interactions  can  be  treated  as  a 
perturbation  (To  61) ,  the  transition  amplitude  is  then  given  by  the  matrix 
element 


T  Q(k 

ocg  o 


V  =  <V")(VIVdb 


(+) 


(k  )> 


2.16 


where  the  and  3  subscripts  represent  the  entrance  and  exit  channels 
respectively.  The  plus  (minus)  superscript  denotes  that,  asymptotically, 
the  incoming  (outgoing)  part  of  the  wave  function  ¥  becomes  identical 
with  the  incoming  (outgoing)  part  of  a  plane  wave.  is  the  inter¬ 

action  potential  between  the  two  nuclei  in  the  exit  channel,  and  U  is 


3 


the  optical  potential  in  the  exit  channel. 


The  DWBA  follows  from  the  observation  that  the  elastic  scattering 
is  usually  the  dominant  process.  This  is  unambiguous  in  the  asymptotic 
region,  but  we  need  the  wave  function  over  all  space  and  particularly  in 
the  interior  region  where  the  projectile  and  the  target  are  strongly  inter¬ 
acting.  The  usual  procedure  is  to  obtain  the  optical  potential  which  repro¬ 
duces  the  observed  elastic  scattering  of  a  on  X  ,  and  which  thus  reproduces 
the  asymptotic  wave  function  accurately.  It  is  then  assumed  that  this  same 
optical  model  wave-functions  may  be  extrapolated  into  the  interior  region  of 
strong  interaction  where  a  and  X  approach  closely.  That  is  we  use  (Sa  65) 

Y  (+)  »  *  X  (+)(k  )  2.17 

a  oc^oc  ac 

and  the  matrix  element  (2.16)  becomes: 

Vk«>y  ■  <yv  ivVw  *  2-18 


where  x’s  are  the  eigenfunctions  of  the  corresponding  optical  potential, 
i.e. 
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{V2+k2  -  ^  [U(r)+U  (r)]}x(k,?)  =  0 


2.19 


where  Uc(r)  is  the  Coulomb  potential. 


Since  the  incident  particle  a  =  b+c  where  b  is  the  outgoing 
particle  and  c  is  the  captured  particle,  then  we  also  have  Y  =  X+c  so 

the  interaction  potential  V  -  U  becomes: 

B  3 


V3~U6  V  "  U&Y  V  +  (V2?X  U£X) 


2.20 


it  is  customary  to  take  V,  as  the  important  interaction  term  for  many 
reasons  (Sa  65,  G1  63,  El  69). 

Consider  a  (d,n)  reaction  then  Eqn.  2.18  can  be  rewritten  as 


T  =  /dr  ,/dr  y  „  ’ (k  ,r  )<Y,n| V__ |x,d>Xj (k. ,r .) 


n  n 


n  n 


np 


'■d  d’  d' 


2.21 


The  central  factor  is  the  matrix  element  of  the  interaction  V  taken 

np 

between  the  internal  states  of  the  colliding  pairs, 

<Y>nl v  Jx,d>  =  fdnip  \jj  V  ip  ip 

np  y  n  np  x  d 


2.22 


it  implies  integration  over  all  internal  coordinates  n  independent  of 

r^  and  '  r^  ,  where  r^  is  the  separation  of  the  centers  of  mass  of  X 

and  d  and  r  i-s  that  Qf  y  and  n  .  This  matrix  element  then  remains 

n 

a  function  of  r,  and  r  .  It  plays  the  role  of  an  effective  interaction 

d  n  r  y 

or  form  factor  for.  the  transition  between  the  elastic  scattering  states 

X(+)  and  X  „  .  It  contains  all  the  information  on  nuclear  structure, 

«  6 

angular  momentum  selection  rules,  and  even  the  type  of  reaction  being  consid¬ 
ered  (stripping,  inelastic  scattering,  knock-on,  etc.).  This  factorization 
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of  the  integrand  of  the  amplitude  2.21  has  the  practical  advantage  that 
one  computer  program  can  treat  a  whole  variety  of  nuclear  models  and 
reaction  modes.  All  we  need  is  a  number  of  optional  subroutines  to 
calculate  the  relevant  form  factors. 


We  can  now  write  the  nuclear  matrix  element  2.22  as 


<y,n  V  X,d>  =  <J  M  ,S  m  V  J„WL,,SJmJ> 
1  np1  yy  nn1  np1  XX  dd 


=  m  <vyV*s  »  (V 

y  y  n  n 


v(r  )  _  (n  n  »r  )i|j  (n  )dn  dn  dnY 

np  rip  ^  P  rip  ^  n  p  X 


2.23 


In  writing  Eqn.  (2.23)  we  assumed  that  the  potential  is  a  central 

potential,  and  the  residual  necleus  Y  is  formed  by  just  adding  a  proton 
to  an  inert  core  (target  nucleus).  The  ru  represents  the  internal 
coordinates  of  particle  i  while  the  r„  represents  the  vector  joining 

the  centers  of  mass  of  i  and  j  .  The  above  assumptions  will  enable  us 

*  -> 


to  expand 


M  ^Vp^PX^ 

y  y 


as 


M  (Vp 
y  y 


.)  =  I<JjMy|jyMy>^(nx)^y(np,?p) 


px 


2.24 


where  <b .  may  be  called  the  captured  proton  wave  function,  and 

jy 

<JiMu|j  M  >  is  the  Clebsch-Gordan  coefficient  as  defined  by  Rose  (Ro  61). 

1  y  y 

n  now  refers  to  the  proton  spin  coordinate.  .  can  be  expanded  further 
p  jy 


as 


WV  “  i£T  "(6p,^)^_ni(np) 


2.25 


because  the  parity  is 


only  those  ^-values  satisfying  Eqn.  2.9  are 


<  -  \  „  —  —  m. 
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allowed.  It  is  customary  to  regard  <f>  as  proportional  to  the  shell 
model  single  particle  wave  function  for  the  orbital  (n£j )  (Sa  65)  so  that, 


R£j (r)  =  ^  unlj (r) 


2.26 


where  u  «.(r)  is  the  normalized  radial  wave  function  and  S„.  is  the 
n  tj  tj 

usual  spectroscopic  factor. 


The  internal  deuteron  wave  function  appears  in  Eqn.  2.23  can 
be  written  as,  (neglecting  the  D-state  of  the  deuteron) 


™  (n  n  ,r  )  =  ip,  (n  n  >r  ) 
S  ,m.  n  p  np  lm,  n  p  np 

d  d  d 


=  <K  (r  )  I  <  -r  m  m  -m  lm  > 
d  np  L .  2  2  nd  n  d 

m 

n 


*P,  (n  )  ip,  (n  ) 

1  i  n  i  i  p 

-o  m'  -7T  m,-m 

2  n  2  d  n 


2.27 


Now  substituting  Eqn.  (2.24)  in  the  nuclear  matrix  element  (2.23)  we  get, 


<JyMy’  ImnlVnplJXMX>1V  =  I<-W|jyMy> 


*  ->  •k 


/<p,(n  ,r  )ip-  (n  )V  (r  )  ip-  (n  n  »r  ) 
J  JM  P  p  l  n  np  np  lm,  n  p5  np 


-an 
2  n 


dnndnp  jf  ^JM(nX^JYM v(nX)dnX 


carrying  out  the  integration  over  nx  we  get, 


:JyV  V 


o  m 
2  n 


V  (r  )  ip-  (n  n  ,r  )dn  dn 
np  np  lm^  n  p  np  n  p 


2.28 


2.29 


Sam  .  I  ^  , _ 


Substituting  Eqns .  (2.25)  and  (2.27)  in  (2.29)  we  get, 


<J  M  ,  i  m  |V  |  JM,  lm  > '  = 
y  y  2  n 1  np 1  XX  d 

l  I  J„M„><:£  b  m  y-m|jy><4  \  m’m  -m'  |  lm  > 

»  xx  yy  2  2  2  n  d  n1  d 


£jmm 


n 


/  m  * 

R«.(r  )[i  Y  P(Q  ,4>  )]  V  (r  )<f>  (r  ) 
^3  P  l  P  P  np  np  d  np 


jh[  (n  )**  (n  )  ^  (n  )  ^  (n  ) 

^U-m  T“n  lK  - - 


2  d  n 


dn  dn 
n  p 


carrying  out  the  integrations  over 
element  reduces  to, 


n  and  n  the  nuclear  matrix 
n  p 


SViMVA’V  = 

l  <JxjV|jyV<4  m  u_m|ju><  1 1  "nVn^d" 

£jm 


9  m  * 

P  £  P  P  np 


where 


m  =  y  +  m  -  m , 
n  d 


and 


D(r  )  =  V  (r  )  cj>,(r  ) 
np  np  np  d  np 


2.30 

2.31 

2.32 


Another  approximation,  concerning  the  potential  V  ,  can  be 
done  here,  namely,  the  zero  range  approximation,  which  consists  of  putting, 
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Now  the  deuteron  wave  function  should  satisfy  the  Schrodinger  equation, 

(v2  -  «2)  *d(s)  =  (M/fi2)D(s) 

and  the  constant  Dq  in  Eqn.  (2.33)  may  be  chosen  to  give  the  correct 
volume  integral  over  D  . 


D 

o 


/D(s)  ds 
2 

f  -  /(V2  -  «2Hd(s)  ds 

cc^2  . 

-  —  /Vs)  ds 


=  -  e/4>d(s)  ds 


where  e 


is  the  binding  energy  of  the  deuteron. 


The  deuteron  wave  function  cj>,  (s)  can  be  taken  as  the 

a 

Hulthen  wave  function. 


<f>d(s)  = 


of  *  q \  -«s  -3s 

<*3  (<*+3)  t  e  -  e 

2tt  (<*-8)  ^ 


Substituting  this  in  Eqn.  2.36  we  get, 


2  4  2  3 

D  I  1.65  x  10  MeV  *F 
o 


The  matrix  element  of  Eqn.  (2.21),  now  becomes: 


T  =  ijS  D°  <JXjVlJyV 


<l\  m  d-m|jP><  3  \  % 

/dr  ;(k  ,r)  u  ».(r)  [i  Y«  (©  *4>)  ] 


rvfj 

(+)  "*■  * 
x'+;(kd,r)  . 


2.34 


2.35 


2.36 


2.37 


2.38 


2.39 
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The  differential  scattering  cross-section  is  then  given  by 

l  |T|2 

do  kn  MAmnmd 

dJJ  =  (2^2)2  kd  (2Jx+1)(2Sd+1) 


2.40 


The  procedure  of  computing  the  transition  matrix  element  of 
Eqn.  2.2l  can  be  summarized  as  follows: 


1)  The  nuclear  matrix  element  is  well  defined  except  for  the  shell 

model  radial  wave  function  u  „.(r)  .  This  is  a  bound  state  wave 

ruij 

function  for  the  captured  particle.  This  wave  function  satisfies  two 
Shrodinger  equations;  the  first  contains  the  nuclear  potential  for  the 
internal  region,  and  the  second, for  the  asymptotic  region, contains  only  the 
coulomb  potential.  The  first  equation  is  solved  with  a  Woods-Saxon  nuclear 
potential  form  with  a  trial  depth  and  this  solution  is  matched  to  the 
solution  of  the  asymptotic  region.  If  the  matching  process  gives  the  correct 
binding  energy  of  the  captured  particle  with  the  target  nucleus,  then  we 
get  the  correct  radial  wave  function.  On  the  other  hand,  if  the  matching 
process  gives  the  wrong  binding  energy.  The  depth  of  the  nuclear  potential 
should  be  changed  and  the  whole  process  is  repeated  until  we  get  the  correct 
binding  energy. 


Once  we  get  the  correct  radial  wave  function  of  the  captured  proton 
the  nuclear  matrix  element  of  Eqn.  2.39  is  completely  defined. 

2)  In  order  to  get  the  distorted  wave  for  both  the  incident  deuteron  and 
outgoing  neutron,  the  corresponding  Schrodinger  equation  of  the  form  given 
in  Eqn.  2.19  should  be  solved.  This  requires  a  knowledge  of  the  appropriate 
optical  potential  parameters.  These  are  usually  fixed  by  fitting  the  elastic 
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scattering  data  of  the  proj ectile- target  system  at  the  same  incident 
energy  and  likewise  the  outgoing  particle-residual  nucleus  system. 

Once  the  optical  potentials  are  defined,  we  can  proceed  to  the 
next  step,  namely,  solving  the  distorted  waves'  Schrodinger  equation. 

3)  To  solve  the  Schrodinger  equations  for  the  distorted  waves,  we  use 
the  usual  partial  wave  expansion  (Sa  64a)  ,  and  the  solution  in  case  of  a 
charged  particle  is  matched  asymptotically  to 

XLJ(k,r)  =  [^(kjr)  -  n^H^^r)]  exp(ipL)  2.41 

where  p  is  the  Coulomb  phase  shift.  Eqn.  2.41  represents  the  solution 
at  large  r  beyond  the  range  of  nuclear  potential  where, 

H^kjr)  =  Gl  +  iFL  2.42 

is  the  outgoing  Coulomb  wave  function  and  is  the  reflection  coeffi¬ 

cient  for  the  (L,J)  wave,  which  can  be  obtained  from  the  matching 
process . 

For  neutral  particles  the  asymptotic  solution  is  the  outgoing 
and  ingoing  spherical  waves. 

4)  Now  after  doing  steps  1,2,  and  3  above  the  transition  matrix  element 
T  can  be  computed,  and  so  is  the  angular  distribution  of  Eqn.  2.40. 

The  above  treatment  is  the  most  simple  treatment  for  computing  the 
angular  distribution  of  a(d,n)  reaction.  Additional  refinements  of  the 
model  have  appeared  in  the  literature,  e.g. 

1)  Use  of  finite  range  potential  instead  of  the  6-function  potential  of 
Eqn.  (2.33), 


>  5  ..  >  ',L*  '  •'  1 
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2)  Use  of  different  forms  of  the  optical  potential  such  as  adding  a  spin 
orbit  term  and/or  surface  or  volume  for  the  real  and  imaginary  parts, 

3)  Inclusion  of  the  d-state  of  the  deuteron,  and 

4)  Use  of  non-local  potentials. 

2 . 3  Compound  Nucleus  Reactions 

As  mentioned  before  (Sec.  2.1)  the  compound  nucleus  model  divides 
the  nuclear  reaction  into  two  steps.  The  first  step  is  the  union  of  the 
colliding  pair  into  a  single  unit,  the  so  called  compound  nucleus .  This 
compound  nucleus,  although  not  stable,  has  many  of  the  properties  of  stable 
nuclei.  In  particular,  it  has  rather  well-defined  energy  levels.  The 
second  part  of  the  reaction  is  the  disintegration  of  the  compound  nucleus, 
either  into  the  pair  from  which  it  was  formed  or  into  another  pair.  This 
disintegration  is  subject  to  probability  laws,  once  the  compound  nucleus  has 
been  formed  the  probability  of  a  particular  mode  of  disintegration  is 
independent  of  the  mode  of  formation. 

The  probability  of  the  formation  of  a  compound  nucleus  is  small 
unless  the  energy  of  the  colliding  pair  coincides  very  closely  with  one  of 
the  energy  levels  of  the  compound  nucleus.  On  the  other  hand,  if  the  coinci¬ 
dence  is  perfect,  the  cross-section  for  the  formation  of  the  compound  nucleus 
shows,  as  a  function  of  the  incident  energy,  sharp  maxima  but  drops  to  small 
values  in  between.  These  sharp  maxima  (called  resonances)  are  ususally  narrow 
and  closely  spaced.  These  resonances  are  individually  of  interest  only  at 
low  energies.  At  higher  energies,  with  increasing  numbers  of  resonances,  the 
physical  interest  (and  usually  the  measurements  itself)  concerns  nuclear  cross- 
sections  averaged  over  the  resonances.  The  statistical  theory  of  nuclear 
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reactions  deals  with  such  average  cross-sections. 

The  treatment  of  average  cross-sections  employs  an  evaporation 
model.  The  evaporation  model  assumes  that  the  compound  nucleus,  in  its 
decay,  loses  all  memory  of  the  way  in  which  it  was  formed,  and  treats  all 
possible  decay  products  as  equally  probable. 

In  the  following  we  will  follow  an  article  by  Vogt  (Vo  68). 


The  evaporation  model  was  first  proposed  by  Bethe  and  Weisskopf, 
in  this  model  a  nuclear  reaction  was  of  the  following  kind: 


Initial  projectile  and  target  nucleus  (labelled  c)  -*  compound 

nucleus  ->  reaction  products  (labelled  c')  is  assumed  to  have  an  average 

cross-section  a  ,  which  may  be  factored, 

cc 


a 

cc 


» 


[a 

comp 


(c  )  ]  [  3  . 


2.42 


The  average  of  the  cross-section  is  over  any  resonance  structure.  The 
first  factor  on  the  r.h.s.  of  Eqn.  (2.42)  depends  only  on  the  initial 
pair  of  particles  c  ,  and  is  called  a  cross-section  for  the  formation  of 
the  compound  nucleus,  the  second  factor  depends  only  on  the  final  reaction 
products  c'  and  is  called  the  branching  ratio.  The  denominator  of  the 
branching  ratio  is  a  sum  of  the  factor  ,  over  all  reaction  products 

(including  c  and  c')  available  to  the  compound  nucleus, 

G  =  I  G  „  2.43 

c»  ° 

The  factorization  assumed  in  Eqn.  (2.42)  is  an  expression  of  the  independence 
of  formation  and  decay  of  the  compound  nucleus. 
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Using  the  receiprocity  theorem  Eqn.  2.42  can  be  written  as 


a  ,  =  a 
cc 


comp  v  ,  2 


k  ,o  (c') 
c  comp 


2.44 


I  kP"aPrmin(cM) 

fi  c  comp 


where  k  is  the  wave  number  for  the  relative  motion  of  the  pair  c  . 

The  evaluation  of  the  denominator  of  Eqn.  2.44  requires  an  enumeration  of 
all  the  possible  reaction  products.  To  carry  out  the  enumeration,  we 
describe  each  reaction  alternative  or  channel  c  (B1  52)  by  the  following 
set  of  quantum  numbers, 


c  =  («,I,i,s,£,J,M  ,tt) 

J 


2.45 


where  «  labels  the  pair  of  particles  in  the  channel  c  represented  by 


"p"  and  their  state  of  excitation  E  where 

P 


E*  =  E  maX  -  E 
P  P  P 


2.46 


E  max  =  E  +  0 
P 


2.47 


where  E  is  the  C.M.  energy  of  the  initial  pair,  E^  is  the  center  of 
mass  energy  of  the  pair  p  ,  and  Q  is  the  Q-value  of  the  reaction.  I 
and  i  are  the  intrinsic  spins  of  the  pair  of  particles  K  ,  s  is  the 
channel  spin, 

")■ 

s  =  I  +  i  2.48 

t  is  the  orbital  angular  momentum  of  the  pair,  J  is  the  total  angular 
momentum, 

J  =  l >  s  =  I  +  I  +  i 

->  -> 

=  J  +  I 


2.49 
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Mj  is  the  z-component  of  J  ,  and  tt  is  the  total  parity.  The  average 
cross-section  a  .  can  now  be  written  as, 

ccoc  f  9 


—  _  _f_  y  (2J+1)  rr  T  / 

.2  /  (21+1)  (2i+l)  T£(  )}< 

K  J  ,  7T 

oc  7 

(«’) 

iTPT }  2-49 

The  unprimed  quantities  refer  to  the  incoming  channels  of  the  reaction, 

\ 

the  primed  quantities  refer*  to  the  outgoing  channels,  and  double  primed 
quantities  in  the  denominator  are  summed  over  all  channels  to  which  the 
compound  nucleus  can  decay.  T^C®)  are  called  transmission  functions 
and  are  given  by  the  optical  model  phase  shifts  6^  of  the  pair  «  , 

2i6  p  2 

T^(oc)  =  1  -  |e  |  2.50 


Those  transmission  functions  are  independent  of  J  and  s  if  the  optical 


model  potential  does  not  contain  spin-dependent  terms,  otherwise  they  should 

be  included  i.e.,  we  replace  each  £  T~(«)  by  £  T  „„(<*)  ,  where 

s  l  1  31  3 
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we  can  now  define  a  channel  wave  function  \jj  ,  describing  all  the  proper¬ 
ties  of  the  channel  except  the  radial  motion, 
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2.52 


where  r 
channel , 
channel , 


is  relative  separation  of  the  pair  of  particles  in  each 


$  is  the  state  of  internal  motion  of  the  two  particles  in  the 

OC 


X  is  the  wave  function  of  the  coupled  intrinsic-spins  I 
sm 

s 
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and  i  of  the  two  particles. 

Outside  the  compound  nucleus  the  whole  wave  function  of  the  system 
can  be  written  as 


¥  ■  I*e*c 

c 


2.53 


where  <J>  is  the  radial  wave  function  in  channel  c  ,  outside  the 
nuclear  potential  the  radial  wave  equation  is 


d  ♦c  ia+1)  ,  2rlckc  2 

—~2 - -  *c  “  — -  *c  +  Vc  ’  ° 

dr  r  c 

c  c 


2.54 


2 

where  n  (=  Z.Z _e  /Tiv  )  is  the  Coulomb  parameter.  The  regular  and 
c  1  z  c 

irregular  solutions  of  Equation  2.45  are  F  and  G  respectively  which 

0  c 

are  normalized  by  their  Wronskian, 


dF  dG 

’  _ £.  _  p  _ 9.  -  v 

'c  dr  c  dr  c 


2.55 


The  outgoing  and  incoming  0  and  I  respectively  can  be  defined  by  a 

C  C 

linear  combination  of  F  and  G  , 


*  ip 

I  =  0  =  (G  -  IF  )  e  C 

c  c  c  c 


2.56 


where  p  is  Coulomb  phase  shift  related  to  the  orbital  angular  momentum 
c 

number  t  and  n  by 
c  c 


P  =  1  tan 

C  n=l 


-i  \ 
n 


2.57 


Asymptotically  I  and  0  have  the  following  behavior, 

C  G 

I*  =  0  =  exp  i  (k  r  -  l  ff-n  ln2k  r  ) 

c  c  r  cc2cc  cc 


2.58 
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The  radial  wave  function  of  Equation  2.54  can  be  written  as  a 

linear  combination  of  I  and  0  , 

c  c 


¥  -  I  —  (A  I  -  BO)* 


c  / v~ 
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c  c  c 


2.59 


where  A  and  B  are  amplitude  coefficients, 
c  c 

I 


corresponds  to  unit  flux. 


v  is  chosen  so  that 
c 


In  terms  of  the  collision  matrix  U  ,  (La  58) ,  the  expansion 

c  c 

coefficients  of  Equation  2.59  are  related  to  each  other  by, 


B 

c 


2.60 


In  a  nuclear  reaction  we  have  an  incident  beam  (plane  wave)  of  a 
given  pair  of  particles  «  ,  and  an  outgoing  wave  of  some  pair  ’  observed 
in  a  direction  (0,<t>)  relative  to  the  incident  beam.  The  total  wave 
function  can  be  written  then  as, 


Y  -  y.  +  ¥  2.61 

me  reac 

The  incident  wave  function  is  described  by 

¥ .  =  $  X  exp(ik  Z  ) 

inc  «  sm  “ 

s 

00  T _  o 

-  $  X  /4tt  7  yim  jo(k  r  )i  y D  2.62 

s  L=0  * 


At  very  large  distances,  when  the  Coulomb  fields  of  the  target  and  projec¬ 
tile  are  completely  screened,  the  radial  part  j^k^r^)  °f  Equation  2.62 

41*  ^ 

can  be  expressed  in  terms  of  the  incoming  and  outgoing  waves  by 
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J<>(k  r  )  =  (k  t  )  1  Fp(k  r  ) 


-  (2v„>  1 


2.63 


Now  equation  2.62  can  be  written  as 
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2.64 


where  the  Clebsch-Gordan  coefficients  <£sOm  |jm  >  and  <^.ls,m«,m  ,  |  Jm  > 

s 1  J  £  s  J 

are  the  coupling  and  decoupling  coefficients  of  the  function  Y^  ^  and 
respectively.  Using  equation  2.52,  Equation  2.64  becomes, 
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2.65 


Since  'J'  __  contains  only  outgoing  waves,  the  coefficients  A  of  the 

1  6aC  Q 

total  wave  function  are  the  coefficients  of  the  incoming  waves  in  Eqn.2.65 


dm  ,  _ _ 

A  =  A  0  =  ik  / irv  <£sOm  I  Jm  >/2£+l 

c  “s t  s'  J 
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for  o'  s  and  m  of  the  incident  beam  and  A  =0  for  all  other  channels 
*  s  c 
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The  amplitude  B  ,  of  the  outgoing  waves  is  given  by  equation  2.60, 
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By  means  of  Eqns .  2.68,  2.65  and  2.61  we  can  write  T  as 
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The  Kroneker-delta  function  6  ,  occurs  in  this  expression  because  the 

incident  plane  wave  contains  outgoing  as  well  as  incoming  waves. 


The  various  possible  cross-sections  are  given  by  ¥ 


reac 


Since  the  actual  detectors  in  the  laboratory  are  presumed  to  select  a 

given  pair  of  particles  in  a  given  state  of  excitation  and  traveling  in 

a  direction  0,<f>,  then  we  can  write  ¥  in  the  form, 
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Comparing  Eqns.  2.69  and  2.68  we  get, 
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s  i  ,  |  Jm  T>  (  5^,  5  ^  ^  ,  6  q  q  ,  —  ^  ,  q  ,  ^  „  q  ) 
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The  differential  cross-sections  for  the  process  «  ->  cc’  is 
given  by  the  following  expression, 
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The  last  expression  can  be  written  in  a  more  convenient  form  (Ba  61) 
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where  PT  (cos  0)  is  Legendre  polynomial  and  B T(oc,cc')  is  given  by, 
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2.72 


2.73 


where  the  Z’s  are  the  Z  coefficients. 
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The  collision  matrix  components  are  complex  numbers  (Vo  68) 
which  vary  rapidly  with  energy.  In  addition  there  is  no  correlation 
between  the  variation  with  energy  of  different  components  of  the  collision 
matrix.  Consequently  when  we  average  the  cross-section  over  an  energy 
interval,  we  can  ignore  all  cross  terms  between  various  components,  so 
each  absolute  square  of  the  collision  matrix  components  takes  the  form: 

uJ 

« ’ s ' V  ,«s£ 
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in  which  we  used  the  optical  model  transmission  functions  of  Equation 
(2.50).  If  we  take  da  ,/dft  from  Equation  (2.72),  making  use  of  Eqns . 

ococ  1 

(2.72)  and  (2.74)  we  get, 
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The  transmission  function  can  be  obtained  from  the  corresponding 
phase  shifts  using  Eqn.  (2.50).  These  phase  shifts  are  obtained  from  the 
solution  of  the  radial  scattering  Schrodinger  Equation  with  the  appropriate 
optical  model  potential.  The  solution  of  these  equations  X ^  (r)  has 
the  boundary  conditions  (Sm  65) 


1)  X£j  (0)  =  0 
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2)  X^(r)  =  a(F£+iG£)  +  b(F£-iG£) 

=  Ft  -  hi  (FriGe} 


with  a  =  j  and  b  =  (1-26^) 


The  transmission  functions  are  then  written  as, 
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As  the  excitation  energy  of  the  compound  nucleus  increases  and 
reaches  the  continuous  region,  or  if  spin,  parity,  and  excitation  energies 
of  the  residual  nucleus  are  not  known,  we  must  include  the  level  density 
expressions  into  any  expression  for  the  transmission  functions.  The 

k 

density  p(E  ,1)  of  states  of  spin  I  of  a  nucleus  is  usually  assumed 

k 

to  depend  on  the  excitation  energy  E  of  the  states  and  their  spin, 
but  to  be  independent  of  the  parity  and  other  properties  (Vo  68) .  It  is 
a  common  practice  to  assume  that  the  energy  dependence  and  spin  dependence 
are  separate  factors  (Vo  68) , 

2 

*  *  -1(1+1) /2a 

p  (E  ,1)  =  p  (E  )[(2I+l)e  ]  2 

where  a  is  called  the  spin  cut  off  parameter,  is  proportional  to 
(E  )T,o  ~  2.5  for  E  =  10  MeV  (Da  68). 
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It  is  found  (Da  68)  that  to  a  good  approximation,  the  p(E  ) 
is  given  by  (Ne  56) , 
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The  any  transmission  function  T^f  should  be  replaced  by 
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where  E^  is  the  maximum  energy  of  the  discrete  region,  E^  =  E 
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CHAPTER  3 

THE  EXPERIMENT 


3.1  The  Experimental  Set-Up: 

The  most  accurate  measurement  of  neutron  energy  is  by  means 
of  time-of-f light  technique.  This  technique  is  characterized  by  good 
resolution  over  a  large  energy  range.  The  idea  is  based  on  measuring 
the  energies  of  neutrons  by  individually  timing  their  flight  over  an 
accurately  measured  path.  It  is  of  great  importance  to  accurately  fix 
two  time  points , 

a)  The  zero-time  point,  which  corresponds  to  the  time  at  which  the 
particle  starts  its  flight  path,  and 

b)  The  stop-time  point,  which  corresponds  to  the  detection  of  the  neutron. 


The  most  general  technique  (Ne  59)  for  obtaining  a  zero  time 
is  that  of  the  pulsed  beam  which  can  give  a  zero  time  independent  of  the 
reaction  mechanism.  In  this  method  the  beam  of  the  incident  particles 
which  initiate  the  reaction  is  pulsed  onto  the  target  for  a  time  short 
compared  to  the  flight  time  to  be  measured  and  at  such  a  rate  that  the 
longest  flight  time  of  interest  can  be  measured  before  the  next  beam 
reaches  the  target.  Since  the  time  of  flight  of  a  several  MeV  neuterons 
is  in  the  nano-second  region,  very  short  ion  bursts  of  high  intensity 
are  required.  Using  both  pre-acceleration  pulsing  and  post-acceleration 
compression,  the  University  of  Alberta  5.5  MeV  Van  de  Graaff  accelerator 
is  able  to  satisfy  these  requirements. 


. 

< 
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An  ion  source  located  within  the  high  voltage  terminal  is  used 
to  ionize  deuterium  gas  by  125  MHz  source  oscillator  and  a  d-c  ion 
beam  is  extracted  by  means  of  a  probe  electric  field.  The  beam  starts 
to  accelerate  downward.  An  Einsel  lens  and  Terminal  Analysis  magnet  are 
used  to  focus  the  beam  which  can  then  be  pulsed  by  two  sets  of  R.F. 
deflection  plates  set  at  right  angles  to  each  other  and  driven  90°  out 
of  phase  by  a  1  MHz  oscillator.  These  R.F.  deflection  plates  sweeps 
the  beam  in  an  elliptical  path  over  a  1/16  inch  aperture,  producing  a 
beam  pulse  of  10  ns  width,  which  after  focusing  by  Focus  3,  is  analysed  by 
a  90°  magnet  to  reflect  the  desired  energy  (momentum).  The  beam  then 
passes  through  a  pair  of  energy  stabilizing  slits  57  inches  away  from  the 
analysing  magnet.  These  slits  feed  the  corona  current  central  system 
which  stabilizes  the  energy.  A  magnetic  quadrupole  lens  is  used  to  focus 
the  beam  inside  the  first  stage  of  the  Mobley  compression  system.  The 
first  stage  of  the  Mobley  system  is  a  pair  of  deflection  plates  fed  with 
a  phase  locked  R.F.  oscillator  of  10  MHz.  Phase  information  is  supplied 
to  the  oscillator  from  a  cylindrical  pick-off  in  the  beam  tube  just  before 
the  R.F.  plates.  The  second  stage  is  a  90°,  75  inch  magnet  whose  focal 
points  are  at  the  target  and  the  R.F.  plates.  The  beam  pulse  is  swept  by 
the  R.F.  field  such  that  the  particles  leaving  the  deflector  plates  first 
take  the  longest  paths  through  the  magnet  while  those  leaving  last  take 
the  shortest  paths.  If  the  phase  and  amplitude  of  the  deflector  are 
properly  adjusted,  all  particles  will  arrive  at  the  target  at  almost  the 
same  time.  A  final  beam  pulse  with  width  of  0.4  ns  can  be  obtained  at 
the  target.  A  schematic  diagram  of  the  beam  handling  system  together  with 
the  Mobley  system  is  shown  in  Fig.  3-1. 
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BEAM  TRANSPORT  SYSTEM 
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In  order  to  determine  the  energy  of  the  outgoing  neutron  groups, 
a  time-of-flight spectrometer  is  used  to  measure  the  neutrons  velocities 
by  timing  the  flight  over  an  accurately  known  distance.  The  stop  signal 
is  obtained  from  a  3  cm.  long  cylinderical  capacitor  placed  30  cm.  before 
the  target.  The  start  signal  is  obtained  from  the  main  neutron  detector, 
which  consists  of  NE  218  liquid  scintillator  in  a  3”  diameter  by  1/2"  long 
quartz  container  optically  coupled  to  an  RCA-4522  photomultiplier  tube. 

The  detector  is  fixed  to  a  remote-controlled  cart.  This  cart  can  be  moved 
through  distances  of  1.0  to  6.3  meters  from  the  target.  The  angle  between 
the  detector  and  the  incident  beam  is  changeable  from  0  to  150°.  In  order 
to  normalize  the  angular  distribution  a  detector  used  as  a  monitor  is  held 
at  a  fixed  position.  The  monitor  detector  consists  of  a  Naton  phosphor 
optically  coupled  to  an  RCA-8575  photomultiplier  tube.  The  monitor  elect¬ 
ronics  system  uses  the  same  stop  pulse  as  the  main  detector  system. 

The  capacitor  pick-off  is  used  to  provide  the  stop  signal  instead 
of  the  detector  pulse,  to  avoid  a  large  number  of  false  starts  (CH  61) 
corresponding  to  beam  bursts  which  do  not  produce  any  neutrons.  This  will 
reduce  the  dead  time  of  the  system. 

Figure  3-2  shows  a  block  diagram  for  the  electronics  of  the 
time-of-flight  spectrometer. 

The  output  of  the  anode  of  the  photomultiplier  tube  is  fed  to  a 
constant-fraction-pulse-height-trigger  (CFPHT) ,  which  triggers  at  1/10  of 
the  input  pulse  height  (Ge  67).  This  was  found  to  eliminate  the  jitter 
(or  walk)  of  the  start  trigger  point  relative  to  the  stop  trigger  point 
due  to  change  in  the  detector  pulse  size.  The  ratio  1/10  is  chosen  for 


FIGURE  3-2 


TIME-OF- FLIGHT  ELECTRONICS 
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the  best  time  resolution.  The  output  of  the  C.F.P.H.T.  provides  the 
start  pulse  for  an  Ortec  model  437  Time-to-Amplitude-Converter  (T.A.C.) 

#1.  The  output  of  the  capacitive  beam  pick  off  is  fed  (after  delayed 
and  amplified),  to  a  Time-Pickof f-Zero-Crossover-Discriminator  (Ge  67a). 

The  output  of  the  Zero-Crossover-Discriminator  is  used  as  the  stop  pulse 
of  T.A.C.  #1  in  the  main  detector  circuit  and  T.A.C.  #3  in  the  monitor 
circuit.  It  is  necessary  for  the  T.A.C.  output  to  be  timed  with  respect 
to  start  pulse,  so  a  signal  corresponding  to  the  start  signal  is  fed  to 
a  gate-and-delay  generator  and  back  into  the  T.A.C.  The  T.A.C.  gives  an 
output  pulse  only  when  the  delayed  start  pulse  is  sensed.  The  output  of 
T.A.C.  #1  is  fed  to  an  Analog-to-Digital-Conver ter  'A'  (A.D.C.A.)  which  in 
turn  feeds  the  S.D.S.-920  on-line  computer. 

The  elimination  of  the  back  ground  y-rays  from  the  neutron 
spectrum,  can  be  achieved  by  measuring  the  time  interval  (using  T.A.C.  #2) 
between  a  delayed  start  pulse  and  cross-over  point  of  an  integrated  and 
clipped  slow  pulse  from  dynode  11  of  the  main  detector.  This  time  interval 
differs  by  about  25  ns  for  neutrons  and  y-rays.  Any  pulse  identified  as 
a  neutron  pulse  derives  a  routing  unit  which  increases  the  address  of  the 
neutron  pulses  to  the  second  2048  channels  while  the  rest  of  the  pulses 
remain  in  the  first  2048  channels. 

For  more  reduction  of  background  the  AiD.C.A.  is  gated  by  a 
coincidence  between  the  delayed  start  pulse  and  the  output  of  S.C.A.  #1. 

The  S.C.A.  //l  is  used  to  set  a  threshold  for  the  main  neutron  detector 


(see  Sec .  3.2). 
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Due  to  the  nature  of  the  monitor  job,  the  high  time  resolution 
and  separation  of  y-rays  are  not  of  great  importance.  So  a  simple  version 
of  main  detector  circuit  is  used  for  the  monitor  circuit.  The  monitor 
spectrum  is  stored  in  a  1024-channel  A.D.C.  and  can  be  dumped  into  the 
computer  memory  on  command.  All  data  are  recorded  on  magnet  tape. 

3.2.  Detector  Efficiency: 

The  detector  efficiency  should  be  known  in  order  to  calculate  the 
absolute  cross-sections.  The  efficiency  of  the  neutron  detector  was  calcu¬ 
lated  using  a  computer  code  written  by  Grandy  (Gr  67). 

The  detector  efficiency  is  a  function  of  the  detector  lower  thres 

hold,  so  the  cutoff  neutron  energy  should  be  measured.  Fig.  3-3  shows  the 

block  diagram  of  the  circuit  used  for  this  purpose.  Fig.  3-4  shows  a 

22 

gamma- ray  spectrum  of  Na  obtained  by  using  the  circuit  of  Fig.  3-3.  The 
lower  threshold  in  electron  energy  is  given  by, 

Lower  Threshold  =  R  •  E 

c 

where  R  is  the  ratio  between  the  cutoff  channel  and  the  channel  corres¬ 
ponding  to  2/3  of  the  height  of  the  0.511  MeV  Compton  edge,  and  Ec=340  KeV 
is  the  maximum  Compton  electron  energy.  To  get  the  threshold  in  terms  of 
proton  recoil  energy  instead  of  electron  energy,  a  response  curve  of  NE  218 
scintillator  to  electrons  and  protons  is  required.  Such  a  curve  is  not 
available  so  the  data  published  by  Smith  et  al.  (Sm  68)  f or  NE  213,  which 
is  quite  similar  to  NE  218,  has  been  used. 
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FIGURE  3-3 

BLOCK  DIAGRAM  OF  CUTOFF  ELECTRONICS 
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FIGURE  3-4 


Na  GAMMA-RAY  SPECTRUM 
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3.3.  Calibration  of  Time  Spectra: 

The  relation  between  the  time-of-f light  and  channel  number  is 
a  linear  one , 


T  =  A  +  B  x  channel  no. 


where  T  is  the  time-of-f light  in  nano-sec.  The  channel  number  of  a 
gamma-ray  peak  along  with  its  flight  time  is  enough  to  determine  A  ,  if 
B  is  known.  The  electronic  circuit  shown  in  Figure  3-5  is  used  to 
determine  B  .  The  T.A.C.  will  give  output  pulses  whenever  the  start- 
stop  time  difference  is  a  multiple  integer  of  the  oscillator  frequency,  and 
the  output  of  the  A.D.C.  will  be  a  chain  of  peaks.  If  d  is  the  average 
spacing  between  any  two  successive  peaks  then  B  is  given  by, 

B  =  lO'Vcd'f)  ns/channel 

where  f  is  the  frequency  of  the  oscillator  in  Mc/sec. 

3.4.  The  ^P(d,n)~^S  Reaction: 

Using  the  experimental  set  up  described  in  Sec. (3.1),  the 
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P(d,n)  S  reaction  has  been  done  at  deuteron  energies  of  4.00  and 
5.45  MeV,  and  flight  distance  of  6.325  m.  A  solid  target  was  prepared 
by  evaporating,  commercially  produced,  Zn^P^  (purity  >  95%)  under  vac¬ 
uum  onto  a  tantalum  backing.  The  monitor  detector  was  placed  at  an  angle 
of  -30°  with  respect  to  beam  direction,  at  3m  from  the  target.  In  each 
experiment  neutron  spectrum  has  been  taken  in  the  angular  range  of  0  to 
140°.  Fig.  3-6  shows  a  typical  neutron  spectrum  taken  at60°  with  4.0 MeV 
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FIGURE  3-5 

BLOCK  DIAGRAM  OF  TIME  CALIBRATION  ELECTRONICS 
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FIGURE  3-6 

NEUTRON  TIME-OF-FLIGHT  SPECTRUM 
31P(d,n)32S  REACTION  AT 
E,  =  4.0  MeV 
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3.5.  Cross-Sections: 


Peak  areas  of  the  5.45  MeV  data  are  extracted  using  a  computer 
code  written  by  J.W.  Tepel  (Te  66).  The  output  of  the  program  is  the  peak 
positions  and  areas  along  with  their  errors.  The  program  fails  to  extract 
the  peaks  of  the  4.00  MeV  data,  this  may  be  due  to  the  low  statistics  of 
this  run,  so  those  peak  areas  are  extracted  by  hand. 


The  absolute  cross-sections  are  calculated  using  computer  code 
written  by  N.  Davison  (Da  69).  This  program  is  written  to  calculate  the 
absolute  cross-section  in  the  C.M.  System  taking  into  consideration  time 
and  monitor  normalization  and  detector  efficiency.  The  program  uses  the 
following  formula  to  calculate  the  absolute  cross-section, 


<§>  =w(i>  . 

abs.  rel. 


mb/sr 


where 


(— ) 
W 


rel. 


m 


W 


1.6  x  A  x 10 
_ m 

Q  n  N 


t .  t . 
1. 1 . 

14 


1 

e 


A  !  A  t  •  t  . 

A  =  A  - — 

m  m  m.£.t. 


where , 

A  =  area  of  a  peak  in  neutron  spectrum 
t.t.  =  true  time  (sec.) 

£.t.  =  main  system  live-time  (sec.) 
e  =  detector  efficiency 
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A'  =  average  of  A' 
m  &  m 

Q  =  average  charge  (y  Coul.) 

ft  =  solid  angle  subtend  by  the  detector 

-4 

=  1.51  x  10  sr  at  6.3  m. 

2 

N  =  number  of  target  nuclei  per  cm 
A^  =  area  of  a  peak  in  monitor  spectrum 
m.£.t.  =  monitor  system  live- time. 


5.8, 


and 


Absolute  cross-section  for  the  G.S.,  2.23,  3.78,  4.29,  5.55, 

32 

6.23  MeV  states  of  S  are  given  in  Appendix  A,  Tables  A-l  to 


A-7 . 
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CHAPTER  4 

ANALYSIS  OF  EXPERIMENTAL  DATA 


4 . 1  Optical  Model  Potentials 


To  carry  out  the  DWBA  calculations  we  need  the  optical  model 
potentials  of  all  particles  involved  in  the  reaction.  These  optical 
potentials  are  expressed  by  a  set  of  parameters.  The  general  form  of 
the  optical  potential  may  be  given  by: 


U(r)  =  V  (r)  -  Vf(x  )  +  i4a  w 
c  v 


.  ..  ,  —  fCx  ) 
w  dx  w 


w 
->  -* 


-(V  +  iw  )  —  ~  f(x  )L*a 

so  so  r  dr  so 


4.1 


V  (r) 
c 


is  the  Coulomb  potential  of  uniformly  charged  sphere  with 


V  (r)  =  ^ 
cV  ;  2R 


[3  -  (~)  ] 
c  c 


r  <  R  =  r  A 
—  c  c 


4.2 


V  (r) 
c  r 


r  >  R 
—  c 


4.3 


The  shapes  of  the  nuclear  potential  wells  are  determined  by  the  Woods- 
Saxon  form  factor 


where , 


f (x)  =  (1  +  eX)  1 


x 

so 


r-r  A 
v 


r-r  A 
so 


so 


x 

w 


r-r  A 
w 


a 

w 


4.4 


4.5 


The  parameters  rv  ,  r^  and  rgQ  are  real,  imaginary  and  spin-orbit 


\ 
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radii  respectively;  the  parameters  a  ,  a  and  a  are  the  real, 
imaginary  and  spin-orbit  diffuseness ' .  The  depths  of  the  potentials  are 
determined  by  the  real  v  ,  the  imaginary  w  and  the  spin-orbit  v 

so 

We  see  from  Equation  4,1  that  it  contains  eleven  parameters. 

These  parameters,  in  general,  are  functions  of  energy  and  atomic  mass 

(Ho  67);  so  they  can  be  determined  by  fitting  the  experimental  elastic 

scattering  angular  distributions  with  those  obtained  from  the  optical 

model  calculations  by  varying  these  parameters  until  a  good  fit  is 

obtained.  Unfortunately,  the  set  of  parameters  which  gives  a  good  fit 

is  not  unique  (Le  64,  Bo  70)  and  this  contributes  to  the  ambiguity  in 

the  deduced  spectroscopic  factors.  So,  to  determine  the  optical  potentials, 

31  32 

necessary  for  the  DWBA  calculations  for  the  reaction  P(d,n)  S  ,  then 
elastic  scattering  data  of  the  different  particles  should  be  available  at 
the  same  energies  and  on  the  same  nuclei.  In  the  following  we  will  discuss 
optical  model  potentials  of  each  particle  individually. 

a)  The  Neutron: 

31  32 

The  reaction  P(d,n)  S  has  been  done  at  E^  =  5.45  and 

4.00  MeV  and  this  yields  neutrons  with  energies  ranging  from  about  12  MeV 

(for  G.S.  at  5.45  MeV)  to  about  4.0  MeV  (for  6.223  MeV  state  at  4.0  MeV). 

32 

There  are  no  elastic  scattering  experiments  for  neutrons  on  S  covering 
this  energy  range.  However,  D.  Willmor  and  P.E.  Hodgson  (Wi  63)  analysed 
a  large  number  of  neutron  elastic  scattering  experiments  on  nuclei  with 
mass  numbers  starting  from  that  of  Si  to  that  of  U  ,  at  energies  from 
1  MeV  to  15  MeV  .  They  deduced  a  general  dependence  of  the  neutron  optical 
potential  parameters  on  mass  number  and  energy.  This  dependence  is  given 
by : 


■ 
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V  =  47.01  -  0.267E  -  0.00118E2  , 

W  =  9.52  -  0.53E  , 

rv  =  1.322  -  0.00076A  +  (4.10_6)A2 

-9  3 

-  (8.10  )A  , 

rw  =  1.266  -  0.00037A  +  (2.10_6)A2  -  (4.10_9)A3  , 

a  =0.66, 
v  * 

a  =  0.48  . 
w 

where  E  is  the  energy  of  the  incident  neutron  and  A  is  the  mass 
number  of  the  target  nucleus. 

b)  The  Deuteron: 

Unfortunately  no  elastic  scattering  data  are  available  for 
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deuteron  with  energies  of  4.0  and  5.45  MeV  on  S  .  So  the  optical 

model  potential  parameters  given  by  D.P.  Gurd  et.  al.  (Gu  68)  have  been 

used.  The  parameters  are  given  in  Table  4-1.  These  parameters  have  been 

28 

extracted  from  elastic  scattering  of  deuterons  on  Si  .  The  set  as 
given  in  Table  4-1  was  used  to  analyse  the  5.45  MeV  data.  A  small  change 
in  w  was  necessary  for  the  4.00  MeV  data  to  keep  the  spectroscopic 
factors  consistent.  The  value  for  w  used  in  this  case  was  25  MeV  . 

c)  The  Proton: 

The  proton  potential  is  a  bound  particle  potential  which  is  a  real 
potential  (see  Sec.  2-2)  .  The  depth  of  the  proton  potential  is  adjusted 
by  the  computer  code  used  for  DWBA  calculations  (see  Secs.  2-2  and  4-3). 


4.6 

4.7 

4.8 
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4.10 

4.11 
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TABLE  4-1 

OPTICAL  MODEL  PARAMETERS 


t 

Parameter 

Deuteron 

* 

Alpha 

* 

Proton 

V(MeV) 

107.0 

200.0 

54.0 

rvCfm) 

1.18 

1.69 

1.25 

a  (fm) 

V 

0.9 

0.576 

0.65 

w(MeV) 

18.5 

15.625 

10.0 

r  (fm) 

w 

1.7 

1.69 

1.25 

a  (fm) 
w 

0.4 

0.576 

0.75 

r  (fm) 
c 

1.25 

1.69 

1.25 

V  (MeV) 
so 

-  - 

-  - 

7.0 

For  Definition  see  Eqs.  4. 1-4. 5 

* 


For  HF  Calculations  Only 
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4 . 2  The  HF  Calculations 

The  contribution  of  the  compound  nucleus  process  to  the  (d,n) 
reaction  has  been  calculated  using  the  statistical  model  by  the  computer 
code  Hauser-Feshbach  written  by  W.R.  Smith  (Sm  65)  and  modified  by 
N.  Davison  (Da  68)  to  include  the  level  density  formula  (see  Sec.  2.3). 

The  output  of  this  program  (differential  cross-section)  was 

very  sensitive  to  any  change  in  the  level  density  parameters  a  and 

b  (eq.  2.80)  of  the  (d ,n)  channel,  while  it  was  completely  insensitive 

to  any  change  in  the  parameters  of  the  other  channels.  The  parameters 

of  the  (d,n)  channel  calculated  by  least  squares  method  and  the  standard 

deviation  in  these  parameters  were  used  to  define  the  region  in  which 

the  HF  calculations  are  quite  certain.  Table  4-2  represents  the  values 

of  the  parameters  of  all  possible  channels  along  with  the  standard 

deviations  a  and  a.  in  a  and  b  respectively, 
a  b 

It  was  also  found  that  the  output  of  the  program  overestimated 
the  magnitude  of  the  differential  cross-section.  To  overcome  this  a 
state  which  shows  a  pure  compound  nucleus  pattern  in  its  angular  distri¬ 
bution  was  used  to  normalize  the  cross-section.  The  normalization  factor 
for  the  4.0  MeV  experiment  was  0.3  and  0.1  for  the  5.45  MeV  experiment. 

4. 3  DWBA  Calculations 

The  DWBA  calculations  have  been  carried  out  by  the  computer 
code  DWUK  written  by  P.D.  Kunz  (Ku  67).  The  optical  potentials 
defined  in  Sec.  (4.1)  with  the  parameters  of  Table  4-1  were  used  to 
carry  out  the  calculations.  The  potentials  are  local  with  a  lower 
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LEVEL  DENSITY  PARAMETERS 
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cutoff  of  zero  fm  in  the  radial  integral  and  a  value  of  0.621  fm  for 
finite  range  corrections. 

Figures  4-1-4-14  show  the  experimental  data  and  both  the  DWBA 

32 

and  HF  predictions  for  the  differential  cross-sections  for  the  S  levels. 

Most  of  the  levels  show  a  strong  direct  pattern  specially  at 
5.45  MeV  deuteron  energy.  The  4.29  MeV  state  is  the  only  state  which 
shows  a  compound  nucleus  pattern  at  both  deuteron  energies.  The  solid 
curve  is  DWBA  normalized  to  the  main  peak.  The  curves  which  are  marked 
with  "X"  represent  the  HF  alone.  The  rhombuses  represent  the  experimental 
cross-section,  and  the  vertical  bars  represent  the  statistical  errors. 

Table  4-3  summarizes  the  DWBA  predictions  for  the  -^-values  of 
32 

S  energy  levels  understudy  along  with  those  obtained  from  previous  (d,n) 

3 

work  carried  out  by  Ferguson  et.  al.  (Fe  68),  (He  ,d)  work  by  Graue  et.  al. 
(Gr  68)  and  Morrison  (Mo  70)  and  those  compiled  by  Endt  and  Van  der  Leun 
(En  67)  before  1967. 

The  figures  of  the  2.237,  4.29  and  5.55  MeV  states  at  4.0  MeV 

deuteron  energy  and  that  of  the  4.29  MeV  state  at  5.45  MeV  deuteron  energy 

have  two  HF  curves.  The  upper  curve  is  the  HF  prediction  with  level 

density  parameters  a- o ^  &nd  b— while  the  lower  one  with  a+o^ 

and  b+o,  . 

b 

4 . 4  Spectroscopic  Factors 

32 

The  spectroscopic  factors  for  S  levels  under  study  have  been 


calculated  according  to  (neglecting  the  interference  between  direct  and 


* 
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P(d,n)  S 

Figure 
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31  32 

P(d,n)J  S 

Figure 

4-3 

31  32 

J1P(d,n)JZS 

Figure 

4-4 

31  32 

P(d,n)  S 

Figure 

4-5 

31  32 

J1P(d,n)JZS 

Figure 

4-6 

31  32 

P(d,nrZS 

Figure 

4-7 

31  32 

P(d,n)  S 

Figure 

4-8 

31  32 

J1P(d,nrZS 

Figure 

4-9 

31  32 

JJ"P  (d  sn)  JZS 

Figure 

4-10 

31_  ,  ,  N 32„ 
P(d,n)  S 

Figure 

4-11 

31-.,  % 32_ 

P(d,n)  S 

Figure 

4-12 

31  {A  x 32 
P(d,n)  S 

Figure 

4-13 

31  32 

J1PCd,n)JZS 

Figure 

4-14 

31^,,  x.  32_ 
P(d,n)  S 

G.S.  angular  distribution  at  E,  =  4.0  MeV 

d 

G.S.  angular  distribution  at  E^  =  5.45  MeV 

2.237  MeV  State  angular  distribution  at  E,  =  4.0  MeV 

d 

2.237  MeV  State  angular  distribution  at  E_  =  5.45  MeV 

d 

3.78  MeV  State  angular  distribution  at  E,  =  4  MeV 

d 

3.78  MeV  State  angular  distribution  at  E^  =  5.45  MeV 

4.29  MeV  State  angular  distribution  at  E,  =  4.0  MeV 

d 

4.29  MeV  State  angular  distribution  at  E,  =  5.45  MeV 

d 

5.55  MeV  State  angular  distribution  at  E^  =  4.0  MeV 

5.55  MeV  State  angular  distribution  at  E^  =  5.45  MeV 

5.8  MeV  State  angular  distribution  at  E,  =  4  MeV 

d 

5.8  MeV  State  angular  distribution  at  E^  =  5.45  MeV 
6.23  MeV  State  angular  distribution  at  E^  =  4.0  MeV 

6.23  MeV  State  angular  distribution  at  E^  =  5.45  MeV 


DIF.  CROSS-SECTION  MB/SR 

0-05  0.10  0.1S  0.20  0.25  0.30  0.35  0.40 


31P  CD,  N)  32S 
ED  =  4.0  MEV 
G  o  S . 

L  =  0  J  =  0 
— *  HF 


30.00  60.00 

C.M.  ANGLE 


90.00 

DEG. 


120.00 


H — 
150.00 


DIF.  CROSS-SECTION  MB/SR 

0.06  0.12  0.18  0.24  0.30  0.36  0.42  0.48 


\ 


DIF.  CROSS-SECTION  MB/SR 

0.04  0.08  0.12  0.16  0.20  0.24  0.28  0.32 


DIF.  CROSS-SECTION  MB/SR 

0.10  0.20  0.30  0.110  0.50  0.60  0.70  0.80 


30.00  60.00 

C.M.  ANGLE 


90.00 

DEG. 


120.00 


H - 

150.00 


DIF.  CROSS-SECTION  MB/SR 

0.15  0.30  O.US  0.60  0.75  0.90  1.05  1.20 


DIF.  CROSS-SECTION  MB/SR 

0.15  0.30  0.45  0.60  0.75  0.90  1.05  1.20 


DIF.  CROSS-SECTION  MB/SR 

O.IO  0.20  0.30  0.40  0.50  0.60  0.70  0.00 


, 


DIF.  CROSS-SECTION  MB/SR 

0.02  0.04  0.06  0.08  0.10  0.12  0.14  0.16 


DIF.  CROSS-SECTION  MB/SR 

0.04  0.08  0.12  0.16  0.20  0.24  0.28  0.32 


DIF.  CROSS-SECTION  MB/SR 

0.08  0.16  0.24  0.32  0.40  0.48  0.56  0.64 


30.00  60.00  90.00 

C.M.  ANGLE  DEG. 


DIF.  CROSS-SECTION  MB/SR 

0.15  0.30  O.HS  0.60  0.75  0.90  1.05  1.20 
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TABLE  4-3 

ANGULAR  MOMENTUM  TRANSFER  VALUES 
FOR  LEVELS  IN  32S  NUCLEUS 


E  MeV 

£ 

P 

X 

Present 

Work 

(Fu  68) 

(Gr  68) 

(Mo  70) 

(En  67) 

G.S. 

0 

0 

0 

0 

0 

2.237 

2 

2 

2 

2 

2 

3.78 

0 

0 

0 

0 

0 

5.55 

2 

2 

2 

2 

2 

6.23 

1 

1 

1 

1 

1 
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compound  mechanisms) , 


(— )  =  (— )  +  (— ) 

W  W 

exp  comp  Dir 


(— ) 


Dir 


■ 


DW 


d  o 

where  (-^-)  is  the  differential  cross-section  given 


DW 


S  is  the  spectroscopic  factor,  and  C  is  the  isospin 


coefficient  which  is  given  by: 


by  DWUK  code, 
Clebsch-Gordan 


C  =  <T  T  T_  T_  T  jT_  ,>  <T_rT  T„TrT_  T_,T_  > 
p  n  Zp  Zn1  d  Zd  X  p  ZX  Zp 1  Y  ZY 


where  T  is  the  isospin,  T  is  the  Z-component  of 

p,n,d,X  and  Y  stands  for  proton,  neutron,  deuteron, 

31  32 

nucleus  respectively.  In  our  case  of  P(d,n)  S  C  ! 


Table  4-4  compares  our  prediction  for  S  and 


T  and  subscripts, 

target  and  residual 
1 

2  ‘ 

those  given  by 


(Gr  68)  ,  (Mo  70) ,  and  (G1  64) . 
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TABLE  4-4 

THE  SPECTROSCOPIC  FACTORS  FOR 
32 S  LEVELS  UNDERSTUDY 


Present  Work 

Experimental 

Theor . 

E  MeV 

X 

V 

4.00  MeV 

5.45  MeV 

(Gr  68) 

(Mo  70) 

(G1  64) 

G.S. 

o+ 

2.63 

2.20 

2.4 

2.2 

2.8 

2.237 

2+ 

1.195 

1.2 

1.3 

1.05 

3.78 

0+ 

0.45 

0.41 

0.4 

0 . 66 

0.47 

5.55 

2+ 

0.17 

0.13 

0.2 

5.78 

1“ 

0.15 

0.11 

0.19 

0.26 

* 

* 

* 

* 

6.23 

U.2) 

0.12 

0.11 

0.12 

0.19 

* 


Calculated  for 


2 


' 
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CHAPTER  5 

DISCUSSIONS  AND  CONCLUSIONS 

It  is  evident  from  Figs.  4-1  -  4-14  that  DWBA  fitting  is 
quite  good  for  all  states  that  show  a  direct  pattern.  The  only  exception 
for  this  is  the  ground  state,  where  the  DWBA  calculations  show  an 
unusually  high  second  peak.  This  may  be  due  to  inadequate  optical 
potential  for  the  deuteron  or  the  outgoing  neutron.  In  spite  of  this 
inadequacy  for  the  ground  state  the  fitting  for  all  other  states  is  quite 
consistent  at  both  deuteron  energies. 

Many  optical  model  parameters,  extracted  from  deuteron  elastic 
scattering  experiments  on  different  neighboring  nuclei  and  in  the  energy 
range  of  4-15  MeV,  have  been  tried  for  calculations.  These  trials  show 
a  large  variations  in  the  deduced  spectroscopic  factors.  This  fact 
stressed  the  need  for  "reliable"  optical  model  parameters. 

It  is  desirable  to  extract  the  optical  model  parameters  from 
deuteron  elastic  scattering  experiment  on  the  same  nucleus  and  at  the 
same  energy  as  those  of  the  (d,n)  experiment. 

Table  4-4  show  that  there  is  a  quite  good  agreement  between  the 
spectroscopic  factors  deduced  from  this  work  and  those  deduced  from 

o 

(He  ,d)  work  (Gr  68,  and  Mo  70)  and  the  theoretical  predictions  of 

28 

Glaudmans  et  al .  (G1  64).  This  agreement  confirms  the  idea  that  Si 

32 

can  be  used  as  a  core  for  describing  the  shell  structure  of  S  . 

The  Glaudmans  et  al.(Gl  64)  shell  model  calculations  has  been 

32 

carried  out  for  several  states  observed  in  S  .  These  calculations  predict 
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the  detailed  wave-function  and  spectroscopic  factors,  bv  assuming  a  Si 
core,  and  employing  two-particle  interaction  between  all  possible  config¬ 
urations  involving  the  extra-core  2S1/2  and  ld3/2  nucleons.  Detailed 

31  32 

wave  functions  of  P  G.S.  and  S  G.S.,  2.237  and  3.78  MeV  states 
are  displayed  in  table  5-1.  Glaudmans  notation  has  been  used  in  this 
table.  This  notation  can  be  explained  as  follows.  The  configuration 
31  1  1 

1031  in  P(J  =  — ,  T  =  — )  ground  state  wave  function  refers  to  two 
2S1/2  particles  coupled  to  J  =  1  and  T  =  0  and  one  ld3/2  particle 
with  J  =  3/2  and  T  =  1/2  ,  so  underlining  indicates  half-integers.  The 
two  terms  are  then  coupled  to  the  final  J  =1/2  and  T  =  1/2.  The 
numerical  entries  are  intensities  intenths  of  a  percent  for  the  various 
pure  wave  functions  at  the  top  of  each  column.  The  sign  refers  to  phases. 


< 
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TABLE  5-1 


Shell-Model  Wave-Functions  Used  to  Calculate 
The  Theoretical  Spectroscopic  Factors 


S3 

S2,dl 

Sld2 

Sld2 

d3 

11 

1031 

1101 

mo 

11 

31P (1/2 , 1/2)  E  =  0.0 

X 

-675 

4 

-288 

-27 

1 

S4 

S2d2 

S2d2 

Sld3 

d4 

00 

0101 

1010 

1111 

00 

32S(0,0)  E  =  0.0 

X 

-496 

391 

59 

-12 

-42 

32S(0,0)  E  =  3.78 

X 

468 

242 

92 

-12 

-186 

S3dl 

S2d2 

S2d2 

S2d2 

Sld3 

1131 

1030 

1010 

0121 

1131 

3 2 S ( 2 , 0 )  E  =  2.237 

X 

I 

-689 

7 

-5 

28 

-240 

Sld3 

d4 

d4 

1151 

20 

20" 

3 2 S ( 2  ,0)  E  =  2.237 

’  X 

8 

20 

-3 
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APPENDIX  A 

TABLES  OF  DIFFERENTIAL 
CROSS-SECTIONS 
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TABLE  A-l 

DIFFERENTIAL  CROSS-SECTION  FOR 
THE  GROUND  STATE 


E,  =  4.00  MeV 
a 

Ed 

=  5.45  MeV 

e° 

W0/dfi)C.M. 

Error 

e° 

(do/dn)C-M. 

Error 

C  .M . 

C  ,M . 

mb/sr 

mb/sr 

mb/sr 

mb/sr 

0.0 

0.351 

0.043 

0.0 

0.836 

0.046 

10.27 

0,270 

0.030 

10.30 

0.402 

0.029 

20.54 

0.079 

0.019 

20.59 

0.067 

0.009 

30.79 

0.152 

0.020 

30.86 

0.096 

0.010 

41.01 

0.180 

0.027 

41.11 

0.070 

0.009 

51.21 

0.126 

0.025 

51.32 

0.209 

0.013 

61.36 

0  o  109 

0.019 

61.50 

0.052 

0.007 

91.57 

0.090 

0.021 

71.62 

0.057 

0.008 

101.55 

0.085 

0.031 

121.50 

0.081 

0.007 

121.36 

0.091 

0.018 

141.11 

0.041 

0.005 

141.01 

0.090 

0.026 
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TABLE  A- 2 

DIFFERENTIAL  CROSS-SECTION  FOR 
THE  2.23  MeV- STATE 


E  =  4.0  MeV 
d 

Ed  ■ 

5.45  MeV 

6° 

C.M. 

(do/d«)c_M_ 

mb/  sr 

Error 

mb/ sr 

e° 

C.M. 

(da/d«)c_M. 

mb/  sr 

Error 

mb/sr 

0.0 

0.291 

0.030 

0.0 

0.581 

0.040 

10.31 

0.262 

0.025 

10.33 

0.420 

0.032 

20.61 

0.433 

0.035 

20.66 

0.486 

0.042 

30.89 

0.275 

0.027 

30.96 

0.592 

0.054 

41.14 

0.316 

0.029 

41.23 

0.463 

0.040 

51.36 

0.242 

0.024 

51.47 

0.342 

0.030 

61.54 

0.235 

0.024 

61.66 

0.118 

0.012 

76.71 

0.205 

0.025 

71.80 

0.218 

0.024 

91.77 

0.251 

0.029 

81.89 

0.160 

0.022 

101.75 

0.178 

0.026 

101.89 

0.210 

0.018 

121.54 

0.248 

0.22 

121.60 

0.149 

0.014 

141.14 

0.233 

0.027 

141.23 

0.182 

0.014 
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TABLE  A- 3 

DIFFERENTIAL  CROSS-SECTION  FOR 
THE  3.78  MeV- STATE 


E ,  =  4.0  MeV 
d 

Ed  ■ 

5.45  MeV 

0° 

(d0/dn>C.M. 

Error 

e° 

(d0/dn)C.M. 

Error 

C.M. 

mb/ sr 

mb/sr 

C.M. 

mb/ sr 

mb/  sr 

0.0 

0.982 

0.076 

! 

0.0 

1.840 

0.092 

10.34 

» 

0.795 

0.061 

10.36 

0.964 

0.060  ! 

i>- 

vO 

• 

o 

CN 

0.381 

0.042 

20.72 

0.316 

0.032  j 

1 

30.98 

0.0  75 

0.023 

41.34 

0.221 

0.026 

41.26 

0.189 

0.030 

51.60 

0.107 

0.014 

51.51 

0.317 

0.035 

61.81 

0.055 

0.014 

61.70 

0.274 

0.029 

71.97 

0.086 

0.018 

76.90 

0.152 

0.027 

121.81 

0.033 

0.002 

91.97 

0.085 

0.021 

141.35 

0.061 

0.016 

101.94 

0.076 

0.026 

121.70 

0.080 

0.019 

141.26 

0.038 

0.024 
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TABLE  A- 4 

DIFFERENTIAL  CROSS-SECTION  FOR 
4.28  MeV-STATE 


E,  =  4.0  MeV 
d 

Ed  - 

5.45  MeV 

e° 

C.M. 

W°/dn;iC.M. 

mb/sr 

Error 

mb/sr 

e° 

C.M. 

(d0/d!2)C.M. 

mb/sr 

Error 

mb/sr 

0.0 

0.148 

0.021 

0.0 

0.077 

0.017 

10.36 

Oo  129 

0.019 

10.38 

0.089 

0,011 

20.70 

0.176 

0.022 

20.74 

0.089 

0.005 

31.02 

0.131 

0.018 

31.08 

0.119 

0.012 

41.32 

0.158 

0.019 

41.39 

0.060 

0.007 

51.57 

0.151 

0.017 

51.66 

0.074 

0.014 

61.77 

0.164 

0.018 

61.87 

0.042 

0.004 

76.98 

0.187 

0.025 

72.03 

0.071 

0.010 

92.05 

0.268 

0.026 

82.13 

0.049 

0.006 

102.02 

0.156 

0.023 

102.13 

0.072 

0.009 

121.77 

0 . 166 

0.017 

121.03 

0.052 

0.008 

141.32 

0.183 

0.023 

141.39 

0.085 

0.008 
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TABLE  A- 5 

DIFFERENTIAL  CROSS-SECTION  FOR 
5  o  55  MeV- STATE 


-  - - - - - ■ -  '  - - -  " 

E  ,  =  4.0  MeV 
d 

— , - , - 

Ed = 

5.45  MeV 

eG 

C  .M  o 

(d0/dfl)C,H. 

Error 

8° 

C.M. 

^d0/dflJC.M. 

i 

Error 

mb/ sr 

mb/ sr 

mb/sr 

mb/sr 

0.0 

0.164 

0.020 

10.41 

0.080 

0.007 

10  „  40 

0.185 

0.020 

20.81 

0.156 

0.015 

20.79 

0.180 

0.022 

31.19 

0.262 

0.033 

31  .  15 

0.190 

0.019 

41.52 

0.219 

0.031 

41 «,  48 

0.214 

0.021 

51.82 

0.185 

0.027 

51.76 

0.124 

0.017 

62.05 

0.083 

0.006 

61.99 

0.142 

0.015 

72.23 

0.082 

0.013 

77.22 

0.097 

0.015 

82.33 

0.054 

0.004 

92.30 

0.118 

0.016 

102.33 

0.048 

0.005 

102.26 

0.105 

0.018 

122.05 

0.071 

0.005 

121.99 

0.131 

0.014 

141.52 

0.085 

0.011 

141.48 

0.155 

0.022 
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TABLE  A- 6 

DIFFERENTIAL  CROSS-SECTION  FOR 
5  .  78  MeV 


E,  =  4.0  MeV 
d 

Ed  = 

5.45  MeV 

0° 

C  .Mo 

(<Wdfi)c_Mi 

mb/sr 

Error 

mb/sr 

6° 

C.M. 

(d0/dQ)C,M. 

mb/sr 

Error 

mb/sr 

0.0 

0.122 

0.018 

0.0 

0.189 

0.023 

10.41 

0.302 

0.025 

10.42 

0.465 

0.018 

20.81 

0.598 

0.041 

20.83 

0.563 

0.026 

31.18 

0 . 460 

0.031 

31.21 

0.476 

0.026 

41.51 

0.453 

0.033 

41.55 

0.217 

0.013 

51.81 

0.234 

0.022 

51.85 

0.249 

0.015 

62.04 

0.266 

0.021 

62.09 

0.154 

0.012 

77.28 

0.198 

0.021 

72.27 

0.197 

0.012 

92.36 

0.161 

0.016 

82.38 

0.095 

0.009 

102.32 

0.144 

0.020 

102.38 

0.043 

0.002 

122.04 

0.209 

0.018 

122.09 

0.119 

0.008 

141.51 

0.183 

0.026 

141.55 

0.135 

0.012 
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TABLE  A- 7 

DIFFERENTIAL  CROSS-SECTION  FOR 
6,226  MeV- STATE 


E,  =  4.0  MeV 
d 

Ea  - 

5.45  MeV 

e° 

C.M. 

(da/dS2)C.M, 

mb/sr 

Error 

mb/sr 

6° 

C.M. 

(d0/dS!)C.M. 

mb/sr 

Error 

mb/sr 

10.41 

0.406 

0.029 

0.0 

0,300 

0.023 

20.81 

0.797 

0.050 

10.44 

0.939 

0.050 

31.18 

0.614 

0.037 

20.86 

0.970 

0.056 

41.51 

0.473 

0.033 

31.25 

0.929 

0.057 

51.81 

0.272 

0.023 

41.61 

0.559 

0.032 

62.04 

0.240 

0.021 

51.92 

0,540 

0.028 

77.28 

0.268 

0.026 

62.17 

0.241 

0.008 

92.36 

0.275 

0.023 

72.36 

0.238 

0.014 

102.32 

0.150 

0.020 

82.47 

0,115 

0,011 

120.04 

0.115 

0.014 

102.47 

0.154 

0.010 

141.51 

0.105 

0.018 

122,17 

0.215 

0.014 

141.61 

0.100 

0,013 
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